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We discuss homogeneous nucleation in a first-order chiral 
phase transition within an efi'ective field theory approach to 
low-energy QCD. Exact decay rates and bubble profiles are 
obtained numerically and compared to analytic results ob- 
tained with the thin-wall approximation. The thin-wall ap- 
proximation overestimates the nucleation rate for any degree 
of supercooling. The time scale for critical thermal fluctua- 
tions is calculated and compared to typical expansion times 
for high-energy hadronic or heavy-ion collisions. We find that 
significant supercooling is possible, and the relevant mecha- 
nism for phase conversion might be that of spinodal decompo- 
sition. Some potential experimental signatures of supercool- 
ing, such as an increase in the correlation length of the scalar 
condensate, are also discussed. 



I. INTRODUCTION 

It is hoped that experiments in the near future at BNL- 
RHIC and CERN-LHC will discover evidence for the so- 
called quark-gluon plasma (QGP) in high-energy colli- 
sions of heavy- ions or protons Q . The QGP is the high- 
density primordial state of strongly interacting matter 
that presumably existed in the early universe a few mi- 
croseconds after the Big Bang pLp). In the QGP, chiral 
symmetry is restored due to finite-teniperature contri- 
butions to the effective potential B (This statement is 
only approximate since finite current-quark masses break 
chiral symmetry explicitly. See below.). 

Given our lack of detailed quantitative knowledge re- 
garding low energy QCD, it is difficult to identify unam- 
biguous observables that would, if found, demonstrate 
beyond reasonable doubt that the QGP had been pro- 
duced. Central ultra-relativistic collisions of heavy ions 
at RHIC and LHC will produce multiplicities on the order 
of 10'' hadrons per event H and raise the possibility of 
analyzing data event-by-event |6|]. Event-by-event fluc- 
tuations might provide an opportunity to discover new 
physics that is usually washed out by standard event av- 
eraging. In particular, if the phase transition is first- 
order with fast expansion so that strong supercooling or 
spinodal decomposition is possible, fluctuations in the 
rapidity density or in the transverse momentum spectra 



of hadrons and isospin fluctuations produced by the for- 
mation of domains of a coherent chiral condensate |0-0] 
could become interesting observables. 

The nature of the QCD phase diagram in the temper- 
ature, T, and baryon chemical potential, /i, plane has 
been intensively studied in recent years ||l^,|ll[ ■ Random 
matrix and effective model calculations suggest that, at 
low /i and non-zero T , a smooth cross-over transition is 
expected and that, at low T and non-zero /i, the chiral 
phase transition should be of first order. Thus, a second- 
order critical point must exist in between these two lim- 
its. The arguments suggesting this topology for the phase 
diagram are in no sense complete, and the phase transi- 
tion could also be first-order for ^ = 0. 

In a first-order phase transition, the thermodynamic 
potential exhibits a metastable minimum corresponding 
to the symmetry restored phase for temperatures slightly 
less than the critical temperature, Tc- This minimum 
gradually disappears as the system cools and ends at 
a point of inflection, i.e. a "spinodal instability". In 
a slowly expanding system, the phase transition would 
proceed through the nucleation of bubbles of the "true 
vacuum" state via thermal activation ||3|,|l^-- [l6| , [l8[ . In 
ultrarelativistic collisions of heavy ions or hadrons, it is 
clear that expansion is very fast compared to that of the 
early universe at the cosmological QCD transition (by a 
factor of ~ lO^^: H ~ lO^^ ^^^-i versus H ~ 10^ sec"!). 
Thus, it is necessary to investigate the time scale for ther- 
mal nucleation relative to that for expansion to see if the 
metastable chirally symmetric state of matter can reach 
the spinodal instability before nucleation is completed. 

In this paper we shall consider a first-order chiral phase 
transition in a rapidly expanding background. We take 
a linear cr-model coupled to quarks as our effective the- 
ory. Treating the quarks as a heat bath which gener- 
ates an effective potential for the soft modes of the chi- 
ral order parameter, we obtain an effective potential at 
the one-loop level. Somewhat below the critical tem- 
perature, this slightly asymmetric effective potential has 
three local minima. Starting with the order parameter 
"localized" in the metastable chirally symmetric state of 
matter, we then study critical bubble profiles and calcu- 
late the nucleation rate to the broken symmetry phase. 
This will provide the temporal and spatial scales of ther- 
mal field fluctuations into the broken symmetry state. 
We will consider both analytical and numerical compu- 



tations. Our numerical results represent "exact" calcu- 
lations. In the analytical treatment, we must turn to the 
thin- wall approximation in order to obtain closed results. 
The comparison between these two approaches will allow 
us to establish limits on the reliability of the thin-wall 
approximation. 

This paper is organized as follows. Section II presents 
our effective theory. Section III describes homogeneous 
nucleation and the methods used for both numerical and 
analytic computations of bubble profiles and nucleation 
rates. Our results are discussed in Section IV. Section 
V presents our conclusions and suggested directions for 
future investigation. We employ natural units, Ti — c — 
ks = 1, throughout. 



mass of the quarks, rUq = gv $^, which enters into the 
expression for the energy, E = ^fW^^^lp^^ . In our phe- 
nomenological approach, we will assume that the quarks 
constitute the heat bath in which the long-wavelength 
modes of the chiral field, i.e. the order parameter, evolve. 
Thus, we will take the effective potential to be 
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To arrive at this form for the effective potential we 
consider a system of quarks and antiquarks in thermo- 
dynamical equilibrium at temperature T (in [ [L1[ Vcff is 
derived also at finite baryon-chemical potential; here we 
shall consider only the case of baryon-symmetric matter). 
The partition function reads 



II. EFFECTIVE THEORY 

As an effective theory of the chiral symmetry break- 
ing dynamics, we assume the linear cr-model coupled to 
quarks i,|Tl,il: 



(2.1) 
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The potential, which exhibits both spontaneously and 
explicitly broken chiral symmetry, is 
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V is the volume of the system. In mean-field approxima- 
tion the chiral fields in the Lagrangian are replaced by 
their expectation values, which we will denote by a and 
7? in the following. Then, up to an overall normalization 
factor. 
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Here q is the constituent-quark field q = {u,d). The 
scalar field a and the pseudoscalar field tt = (711,772,773) 
together form a chiral field <& — {a, tt) . The param- 
eters of the Lagrangian are chosen such that chiral 
SUl(2) (E) SUii{2) symmetry is spontaneously broken in 
the vacuum. The vacuum expectation values of the con- 
densates are {a) = /„ and (7?) = 0, where /^ = 93 MeV is 
the pion decay constant. The explicit symmetry breaking 
term is due to the finite current-quark masses and is de- 
termined by the PC AC relation which gives hg = /ttTiJ, 
where m.^ = 138 MeV is the pion mass. This leads to 
^2 = /2 _ j„2/;y2_ rpj^g yj^i^g Qf ;^2 ^ 20 Icads to a a- 



2A^/^ 



equal to 600 MeV. In mean 



mass, m^ 

field theory, the purely bosonic part of this Lagrangian 
exhibits a second-order phase transition [£1| at Tc = \/2u 
if the explicit symmetry breaking term, hq, is dropped. 
For hq ^ 0, the transition becomes a smooth crossover 
from the restored to broken symmetry phases. For g > 0, 
the finite-temperature one-loop effective potential also in- 
cludes the following contribution from the quarks: 
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Here, dq — 24 denotes the color-spin-isospin-baryon 
charge degeneracy of the quarks. We note that Vq{'^) 
depends on the order parameter $ through the effective 
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where 
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Taking the logarithm of Z, the determinant of the Dirac 
operator can be evaluated in the standard fashion [ p2| , 
and we finally obtain the grand canonical potential 
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where 
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For our purposes, the zero-temperature contribution to 
Vq can be partly absorbed into U via a standard renor- 
malization of the bare parameters A^ and u^. However, a 
logarithmic correction from the renormalization scale re- 
mains and is neclected in the following (see [|ll| for more 
details). The finite-T part defines Vq as given in eq. (p!^). 
We now turn to a discussion of the shape of the effec- 
tive potential at various temperatures. 




FIG. 1. The one-loop finite-T effective potential for 



Tsp < T < Tc and coupling constants g = 5.5, A 
The curves are labeled by the temperature in MeV. 
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For sufficiently small g one still finds the above- 
mentioned smooth transition between the two phases. At 
larger coupling to the quarks, however, the effective po- 
tential exhibits a first-order phase transition. For tem- 
peratures near the critical temperature, T4ff displays a 
local minimum $ = ^f{T) ~ which is separated by a 
barrier from another local minimum at $ = $t(r) > 0. 
(There is another local minimum for negative $ which is 
of higher energy and need not concern us.) These two 
minima are degenerate at T = Tc. For example, g = 5.5 
leads to a critical temperature of Tc = 123.7 MeV. Unless 
stated differently, throughout the manuscript we shall 
adopt the values g ^ 5.5 for the quark- field coupling and 
A^ = 20 for the self-coupling of the chiral fields. 

Chiral symmetry is (approximately) restored for T > 
Tc. The minimum at $ = <i>t becomes the absolute min- 
imum as the temperature is reduced below Tc. As the 
temperature is lowered, the local minimum at $ ~ ap- 
proaches the intervening maximum. These two extrema 
meet and form an inflection point at the spinodal tem- 
perature, Tsp. Only the broken symmetry minimum re- 
mains for T < Tsp. For our standard choice of parame- 
ters, Tsp = 108 MeV. The potential in the cr-direction is 
shown in Fig.Jl for the range of interest, Tgp < T < Tc; 
see also Ref. M\ Note the rather small barrier between 
the two local minima of Ves at Tc. The first-order transi- 
tion for g = 5.5 is weak compared to those conventional 
in bag model calculations |lq|. This weakness is due to 
the fact that gluons, absent in the present linear cr-model, 
are treated as free fields (for T > Tc) in the bag model. 
It is not our goal here to argue which approximation is 
more reasonable. Rather, we view our approach as com- 
plementary to previous studies which employed the bag 
model [^|12|jl6f|. The existence of a local minimum in 



^The right diagram in Fig. 1 of Ref. M is incorrectly indi- 
cated as corresponding to T = 100 MeV. It actually represents 
T = 108 MeV. 



Vcs raises the possibility of super-cooling effects since the 
true broken symmetry minimum can only be reached by 
tunnelling. The weakness of our first-order transition is 
likely to overestimate tunnelling rates and underestimate 
the effects of supercooling. 
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FIG. 2. The one-loop finite-T effective potential for 
Tp < T < Tc and coupling constants g = 10, A^ = 20. The 
curves are labeled by the temperature in MeV. 

In fact, increasing the coupling to the quarks leads to 
a larger barrier between the two local minima of the ef- 
fective potential, as seen in Fig. |[ In other words, the 
first-order transition becomes stronger, with larger la- 
tent heat. Furthermore, the spinodal temperature Tsp 
decreases considerably, to about ~ 78 MeV. As will be- 
come clear below, both effects (the higher barrier and 
the lower temperatures) act to reduce the decay rate of 
the metastable minimum as compared to our "standard" 
choice of couplings corresponding to the potential shown 
in Fig. |I|. 
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FIG. 3. The one-loop finite-T effective potential for 
Tsp < T < Tc and coupling constants g = 5.5, A'^ = 2.2. 
The curves are labeled by the temperature in MeV. 

There is also the possibility that the cr-meson is much 
lighter at Tc than at T = Q (i.e. A^ is small) if QCD is 
near a chiral critical point. Chosing for example A^ ~ 2.2 



with the pion decay constant /^ and vacuum mass m^r 
fixed, such that v'^ = 0, leads to the potential shown in 
Fig. ^. Clearly, decreasing A^ in this theory reduces Tc 
very much, to less than 100 MeV. Reducing g has practi- 
cally no effect, except that the first-order phase transition 
disappears at g < 4. Such low Tc appear to be excluded 
by present lattice QCD results |17|, and will therefore 
not be considered further. 

In order to obtain approximate analytic formulas, it is 
convenient to express Vcs over the range < $ < T in 
the familiar Landau-Ginzburg form (see also p^) 
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This form is adopted purely for reasons of convenience 
and no deeper meaning should be attached to it. It is, 
for example, obviously incapable of reproducing all three 
minima of VeS- We shall use this expansion both for a 
limited range of temperatures (Tgp < T < Tc) and for 
a limited range of field strengths (0 < <& < T) which 
includes the two local minima of interest. With these 
restrictions, this polynomial form is found to provide a 
quantitative description of Vcff- The coefficients a„ are 
determined for each T by performing a least squares fit 
to the true effective potential. 



III. HOMOGENEOUS NUCLEATION 



We wish to consider the physical situation (realizable 
for Tgp < T < Tc) in which a scalar field, corresponding 
to the (T-direction, is initially localized in the metastable 
restored symmetry minimum of an asymmetric double- 
well potential. (We neglect the higher-energy local min- 
imum with <I> < 0.) It is well known that such a system 
can decay by the nucleation of thermally activated bub- 
bles of the true vacuum inside the false one P,p"^- p^ , p^ . 
The nucleation rate per unit volume per unit time is ex- 
pressed as 



-F,/T 



(3.1) 



where -FJ, is the free energy of a critical bubble and where 
the prefactor V provides a measure of the saddle point 
of the Euclidean action in functional space. It is conven- 
tional to write P as a product of the bubble's growth rate 
and a factor proportional to the ratio of the determinant 
of the fiuctuation operator around the bubble configura- 
tion to that around the homogeneous metastable state. 
The following analysis will concentrate on the exponen- 
tial barrier penetration factor, and we will approximate 
V hy T^. This approximation is known to give an overes- 
timate of r, see e.g. |l3,M. The reason for this focus is 
that Fh is necessarily infinite at Tc and zero at Tgp. The 
exponential factor is thus expected to make the dominant 
contribution to the structure of T. 



In order to determine the role of bubble nucleation in 
the evolving system and to be able to compute the de- 
cay rate r, it is necessary to study the critical bubble 
and some of its features. The critical bubble is a radi- 
ally symmetric, static solution of the Euler-Lagrange field 
equations that satisfies the boundary condition <I>(r — > 
oo) — > <!>/ where $/ is the false vacuum of the effec- 
tive potential. Energetically, this boundary condition 
corresponds to an exact balance between volume and 
surface contributions which defines the critical radius, 
R = Re- The critical bubble is unstable with respect 
to small changes of its radius. For R < Re, the surface 
energy dominates, and the bubble shrinks into the false 
vacuum. For R > Re, the volume energy dominates, and 
the bubble grows driving the decay process. 

The critical bubble can be found by minimizing the 
free energy 
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r dr 
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with respect to the field <I>. This leads to 
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where 
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(It has been assumed in Eq. (3.2) that the effective po- 
tential is shifted so that Vcs{^f,T) = 0.) The desired 
solution of Eq. (3.3) is found by imposing the boundary 



condition <f>(r -^ cx)) > $/. 

For numerical purposes it is more convenient to impose 
boundary conditions at r = by setting <I>'(0) = and 
$(0) equal to an "escape value" $e- It is straightforward 
to integrate Eq. (3.3) to obtain the bubble profile <i>(r) 
for a given temperature. The escape value is then ad- 
justed until the desired boundary conditions at r -^ cx3 
are realized. This process can be facilitated by smoothly 
joining the numerical values of <^(r) to the analytic so- 
lution to Eq. ( |3.3D that is regular as r ^ oo. Close to 
the critical temperature, $(?■) has a simple qualitative 
behaviour. It stays close the escape value $e for r < Re 
and then makes a rapid transition to the false vacuum, 
$/. For obvious reasons. Re is identified as the radius 
of the critical bubble. Once the profile of the critical 
bubble has been found, the free energy Ff, is obtained by 
computing the integral in Eq. (3.2) numerically. Finally, 
the decay rate can be calculated using the approximate 
prefactor V = T*. 

A simple physical analogue of this procedure can be 
obtained by turning the potentials in Fig. [| upside down 



and considering the mechanical motion of a point par- 
ticle moving in the inverted potential. The problem of 
determining the critical bubble profile is then seen to be 
equivalent to finding that point on the true vacuum hill 
such that a ball start ing at rest and obeying the equa- 
tion of motion, Eq. (|3.3| ), would descend and land at 
rest on the top of the false vacuum hill at the point <&/ . 
The starting point is evidently the escape point. It does 
not correspond to the true vacuum point but rather to 
some smaller value. This mechanical analogy also makes 
it clear that —Vcs{^e) must be greater than —Vcs{^f) 
due to the presence of the second, "frictional" term in 
Eq. (p^). As we shall see, the critical radius Re di- 
verges as T -^ Tc- The frictional term then becomes 
small in some sense. The various approximate methods 
for dealing with this frictional term approximately are all 
described as "the thin- wall approximation" . 

Although the thin-wall approximation is not expected 
to be valid near the spinodal, it should provide a reliable 
description of nucleation near Tc- Given the frequency of 
its adoption and its simplicity, we find it useful to apply 
the thin-wall approximation to th e present problem. A 
quartic potential such as Eq. ( 2.10D can always be rewrit- 
ten in the form 



W{(p) ^a {ip^ - c?f + jif 
Specifically, we have 
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The new potential W{(p) reproduces the original Vcs{^) 
up to a shift in the zero of energy. We are interested in 
the effective potential only between Tc and Tgp. At T^ we 
will have two distinct minima of equal dept h. This clearly 
corresponds to the choice j = in Eq. (3.5) so that W has 
minima at ip ^ ±a and a maximum at (^ = 0. The mini- 
mum at if — —a and the maximum move closer together 
as the temperature is lowered and merge at Tgp- Thus, 
the spinodal requires j/aa^ = — 8/3\/3 in Eq. (3^). In 



the present case, the parameter a is essentially indepen- 
dent of temperature, and a rises roughly linearly by some 
3% as the temperature falls from Tc to Tgp. The parame- 
ter j/aa^ falls roughly linearly from (at Tc) to — 8/3V3 
(at Tsp). 

The explicit form of the critical bubble in the thin- wall 
limit is then given by |Tq] 



Vc 



■,^,Rc) = (Pf 



1 



^V2^ 



1 — tanh 



Re 



where (pf is the new false vacuum, Re is the radius of 
the critical bubble, and £, = 2/to with m? = W"{(pf) 
is a measure of the wall thickness. The thin- wall limit 
corresponds to ^/Rc ^ 1 [@, which can be rewritten 
as (3|j|/8aa'^) <C 1. This small parameter has the value 
of l/VS at the spinodal, which suggests that the thin- 
wall approximation might be qualitatively reliable for our 
purposes. Nevertheless, it merits a quantitative check. 
In terms of the parameters a, a, and j defined above, we 
find 
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in the thin- wall limit. Determination of the critical radius 
requires the surface tension, E, defined as 

The critical radius then becomes 

2S 
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(3.5) where 



AW = V{^f)-V{^t)^MJ\ ■ (3-15) 

The energy of a critical bubble is finally given by 
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(3.16) 



From knowledge of Fi, = Ec, one can evaluate the nucle- 
ation rate F. In calculating thin-wall properties, we shall 
use the approximate forms for (j>t, (f>f, S, and AW for all 
values of the potential parameters. 



IV. RESULTS 

We begin by showing the critical bubble profiles ob- 
tained from the exact numerical solution of Eq. (B^) in 
Fig. |. 
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FIG. 4. Exact bubble profiles for temperatures 110, 112, 
114, 116, 118, 120, and 122 MeV. Bubble size increases mono- 
tonically with temperature. Recall that Tc — 123.7 MeV and 



T,p = 108 MeV. 



For small supercooling the bubble profiles exhibit a 
pronounced "core" with $ ~ $( separated by a rela- 
tively thin wall from the region $ = $^. On the other 
hand, critical bubbles become "coreless" for stronger su- 
percooling; their thickness is of the same order as the 
radius, and the field at the escape point, $(r — 0), does 
not reach <I>t. For comparison. Fig. |5| shows bubble pro- 
files for temperatures T = 110 and 12 M eV as obtained 
from the numerical solution of Eq. ( p.3[ ) and from the 
thin-wall approximation. In this case, we have used the 
same polynomial approximation to VcS in both exact and 
thin- wall calculations so that discrepancies can be associ- 
ated cleanly with the thin- wall approximation. (The crit- 
ical bubble profile is quite sensitive to the precise form of 
Vef[. This is particularly true when T Ki Tc and the bub- 
ble radius is large.) As expected, the thin- wall approxi- 
mation provides a qualitatively reasonable description of 
the critical bubble profile. 




FIG. 5. Critical bubble profiles for temperatures 110 and 
120 MeV. The solid line indicates exact numerical results. 
The dashed line shows the results obtained with the thin-wall 
approximation. 
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FIG. 6. Critical radii computed numerically (solid line) and 
with the thin-wall approximation (dashed line) as functions 
of the degree of supercooling. 

In Fig. g, we compare the critical radii Re obtained 
from the two calculations as a function of temperature. 
For the exact numerical calculation. Re has been ex- 
tracted from the bubble profile shown in Fig. ^ as the 
oint at which the curvature of the profile changes sign 
This figure makes it clear that the thin-wall approx- 
imation becomes quantitatively correct as T —> Tc and 
i?c — > cxD as expected. The fact that the thin-wall ap- 
proximation systematically underestimates Re will be re- 
flected in a systematic overestimate of the decay rate of 
the metastable state. 

This expectation is also supported by Fig. H which 
shows Fh/T as a function of the degree of supercooling, x, 
defined as x = {Te — T)/{Te — Tsp). This figure also allows 
us to calibrate the temperature range of primary interest 
in the present work. We have chosen to concentrate on 
the relativ ely v iolent behaviour of the exponential fac- 
tor in Eq. (3.1) and to treat the prefactor, V, crudely. 
Given, the dramatic changes in F^/T as a function of 
supercooling shown in this figure, this focus seems ap- 
propriate. For temperatures such that Fb/T > 1, F will 
be strongly suppressed by the exponential factor, and 
the system is likely to remain in the metastable state for 
times comparable to the expansion time of the system. 
For lower temperatures, the exponential factor is unim- 
portant, and r is determined primarily by the prefactor, 
V . At such temperatures, it becomes necessary to make 
a more detailed comparison of the decay rate with the ex- 
pansion rate of the system in order to decide the fate of 
the metastable state. Such comparison requires a more 
reliable description of V than that adopted here. The 
exact results shown in the figure suggest that the system 
is likely to supercool for some two-thirds of the way from 
Te to Tgp- By contrast, the thin- wall results yield the 
significantly smaller value of a; w 0.4. 



This definition of Re is somewhat arbitrary. This figure 
is intended only for the comparison of exact and thin-wall 
results. Our numerical calculations do not employ Re in any 
of the quantities discussed below. 
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FIG. 7. The free energy of the critical bubbles as a func- 
tion of the degree of supercooling. The solid line shows the 
exact (numerical) result, the dashed line shows the thin- wall 
approximation, and the dotted line shows exact results for the 
modified effective potential discussed in the text. 

Given the weakness of the first-order phase transition 
predicted by our model with g ~ 5.5, it is remarkable that 
the system can come so close to the spinodal temperature 
before the stability of the metastable state is seriously 
challenged. It is to be expected that effective potentials 
which maintain the same Tc and Tgp but which have a 
higher barrier between the phases at Tc would permit 
the metastable state to come even closer to the spinodal 
temperature. 
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FIG. 8. The effective potential at T = 116 MeV. The solid 
line shows the original 1/cft from Fig. gj. The dotted line shows 
the modified VcH as discussed in the text. 

This point can be illustrated by considering an artifi- 
cially modified effective potential with a higher barrier. 
We construct this potential from our quartic approxima- 
tions by arbitrarily multiplying the parameter j at each 
temperature by a factor of {T^ — T)/{Tc — Tap)- Evi- 
dently, this modification has no effect on Ves at the crit- 
ical temperature (where j is 0) or at the spinodal (where 
the additional factor is 1). The maximum difference in 
effective potentials lies halfway between these extremes 
and is illustrated in Fig. K While the difference between 



these potentials is not dramatically large, it is sufficient 
to raise Fb/T from the value of 2 (for the original Ves) to 
10. Following exponentiation, this leads to a reduction 
of r by a factor of almost 3000. The values of Fb/T for 
this modified potential, shown in Fig. M, suggest stability 
of the metastable phase for all x < 0.83. 

It is now a simple matter to obtain F, the nucleation 
probability per unit time and volume, and to determine 
the space-time scales for homogeneous nucleation. This 
decay rate must be compared with the expansion rate of 
the central region of a high-energy collision. Certain ob- 
vious conditions apply. For example, the critical radius 
Re cannot exceed the horizon. Otherwise, the center of 
the bubble cannot be causally connected to its surface, 
and the hydrostatic balance of surface and volume energy 
is thus impossible. In addition, there is a time scale asso- 
ciated with a spontaneous coherent thermal fluctuation 
of energy Ec with a length-scale of Re- If that time-scale 
is larger than the inverse expansion rate, bubble nucle- 
ation can be regarded as a slow process, and deep su- 
percooling is to be expected. (This is the condition that 
Fh/T ^ 1 discussed above.) To our knowledge, these 
questions have not previously been addressed in studies 
of homogeneous nucleation in the QCD phase transition 
as found in hadronic collisions. 

To simplify our considerations, we assume that the cen- 
tral rapidity region exhibits one-dimensional (longitudi- 
nal) Hubble flow ||2J| with an expansion rate (i.e. Hubble 
constant) H = 1/t, where t denotes proper time. Even in 
large nuclei, transverse expansion can increase the local 
expansion rate considerably when the system has cooled 
down to r w Tc p5||. Reasonable expansion rates are 
H fn 0.1-1 fm-i ^^. 
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FIG. 9. r^'* divided by an assumed expansion rate of 
H = 1 fm~^ as a function of temperature. The solid line 
shows the exact (numerical) result, the dashed line shows the 
thin-wall approximation, and the dotted line shows exact re- 
sults for the modified potential. 

We have chosen r^''^/H as a measure for the compar- 
ison of nucleation and expansion |13| . The results shown 
in Fig. Q assume an expansion of H = 1 fm~^. When the 



nucleation rate is much larger than the expansion rate, 
it is reasonable to adopt standard thermodynamic tech- 
niques and to describe the phase transition by an ideal- 
ized Maxwell construction [|l6[. Our results suggest that 
the Maxwell construction is not likely to provide a reli- 
able description of the phase transition dynamics. While 
the results of Fig. do not provide an unambiguous pre- 
diction for the fate of the metastable state, they do make 
it appear likely that the rapidly expanding system has 
an appreciable probability of remaining in the restored 
symmetry phase even close to the spinodal instability. 
Thus, at least some fraction of all heavy ion events should 
show traces of this non-equilibrium transition. This un- 
derscores the importance of the event-by-event analysis 
of data g . 

We now consider some potential consequences of a non- 
equilibrium transition associated with strong supercool- 
ing relative to an adiabatic transition. First, the charac- 
ter of the phase transition will certainly influence global 
observables and the expansion of the hot system. For 
example, a non-equilibrium transition is always associ- 
ated with entropy production which can be as large as 
30% |1^. This would be reflected in the multiplicity 
of hadrons in the final state. Moreover, phase coexis- 
tence in equilibrium leads to a distinct hydrodynamic 
expansion pattern which would be absent in a rapid non- 
equilibrium transition with strong supercooling. Since 
this issue has already received considerable attention in 
the literature p6| , we refrain from a detailed discussion 
here. 

Further, as discussed in |^, if the chiral condensate is 
"trapped" in the false minimum and supercooling per- 
sists close to the spinodal instability, the subsequent re- 
alignment in the chiral field in the direction of the true 
vacuum can generate strong coherent pseudoscalar fields, 
7?, which will eventually decay into a "beam" of coherent 
pions. The experimental situation will then be similar to 
that following an "instantaneous quench" of the restored 
symmetry phase ||8|. 

Consider the behavior of the correlation length of the 
^-field. It is obvious from Fig. |l| that the curvatures of 
the effective potential at the minima $/ and <I>t are rather 
different for T < Tc- Thus, the effective mass, m^jj — 
d^Vcff/d$^, and the correlation length, I — l/rriofr, of 
the field will also be different. In particular, if the rapidly 
expanding system supercools appreciably and if the order 
parameter is trapped in the metastable state, one can 
expect a clear increase of I. 
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FIG. 10. Correlation length of the chiral order parameter 
in the metastable state (If) and in the global minimum corre- 
sponding to an equilibrium transition (It). Solid lines corre- 
spond to the exact results, while dashed lines were obtained 
within the thin- wall approximation. 

In Fig. |l^ we show I at $/ and at $t as a function 
of the degree of supercooling. The correlation length at 
<I>j increases as the metastable minimum disappears. By 
contrast, there is a smooth decrease in the correlation 
length at the true minimum $( (which would be the ex- 
pectation value of $ if the transition proceeded in equi- 
librium). In the present model If can exceed It by as 
much as a factor of 2-3, depending on the degree of su- 
percooling. The precise values of If and It depend on 
the specific effective model adopted, but the qualitative 
observation that If increases dramatically as the system 
approaches the spinodal instability is general. 

At the critical temperature, the thin-wall approxima- 
tion leads to only one correlation length, Ztw = C/2, where 
^ is the thickness of the critical bubble as given by Eq. 
(3.12). This is because the thin-wall approximation as- 
sumes both a small degree of supercooling and a small 
latent heat, i.e. a weak transition. Thus, the curvatures 
at $/ and at $t coincide to leading order in £,/Rc- To 
leading order in (,/ Re and for Tgp < T < Tc, we find 



W"{^f)-W"i^t) = ~6\j\/a 



(4.1) 



Correlation lengths in the thin-wall approximation are 
shown in Fig. 11^ by the dashed lines. The splitting of If 
and It is described reasonably well (at least qualitatively) 
except close to the spinodal, where this approximation is 
most unreliable. 

We note that the horizon, 1/H, in expanding systems 
provides an upper bound on the correlation length; If 
cannot increase beyond the horizon. Thus, observable 
consequences can be expected only if the expansion rate 
diminishes near the spinodal. On the other hand, corre- 
lation lengths cannot become large if expansion rates are 
too small, e.g. H < 0.1 fm~^, since the dynamics will no 
longer support strong supercooling. 

As we have argued above, nucleation is not effective for 
small supercooling because Fb/T is large. The decay rate 
is then likely to be much smaller than the expansion time. 



and the system is likely to remain in the false vacuum. 
On the other hand, the approximations of homogeneous 
nucleation theory are likely to break down near the spin- 
odal where the decay rate of the false vacuum becomes 
large. There, the condition R,. ^ Zy, inherent to nucle- 
ation theory, is no longer valid. If Re ~ If, it makes 
little sense to think of a bubble as a coherent superposi- 
tion of random thermal fluctuations of wavelength If in 
hydrostatic equilibrium with the surrounding symmetric 
phase. It may rather be more appropriate to treat fluc- 
tuations on the scale If as random thermal noise to be 
introduced into the classical equations of motion for the 
order parameter through, for example, a Langevin equa- 
tion 1^. In short, our results suggest that homogeneous 
nucleation theory may never provide a suitable modelling 
of the chiral phase transition in particle collisions. 



V. SUMMARY AND OUTLOOK 

We have considered homogeneous bubble nucleation 
in a first-order chiral phase transition within an effective 
field theory approach to low-energy QCD, i.e. a linear 
cr-model coupled to quarks. Integrating out the quarks 
leads to an effective potential for the soft modes of the 
chiral fields (i.e. the order parameter) with two compet- 
ing minima for T = Tc, corresponding to a first-order 
phase transition. The approximation of this effective po- 
tential by a familiar Landau-Ginzburg form allowed us to 
obtain analytical results for critical bubble profiles and 
decay rates using the thin-wall approximation. Compar- 
isons with exact numerical solutions, which do not em- 
play the thin-wall approximation, suggest that the thin- 
wall approximation is not adequate for the description of 
high-energy hadronic collisions. Its use results in a signif- 
icant underestimation of the size of critical bubbles. The 
effects of this error are then amplified in the calculation 
of the decay rate. Furthermore, the thin-wall approxi- 
mation can not be used to obtain the correlation length 
in the metastable state close to the spinodal. 

Our results show a rapid decrease with temperature 
in the free energy of critical bubbles from infinity at the 
critical temperature to zero at the spinodal. Given the 
exponential dependence of T on Fh/T, this fact provides 
some justification for our crude treatment of the prefactor 
in r. Furthermore, it suggests that the metastable phase 
is likely to survive at temperatures such that Ft ^ T 
is independent of the details of either the decay pro- 
cess or the expansion. Even the weak first-order tran- 
sition considered here suggests the viability of the re- 
stored symmetry phase for (T^ - T)/{Tc ~ Tsp) < 2/3. 
As we have suggested, stronger transitions will allow the 
metastable state to survive even closer to the spinodal in- 
stability. The fate of the metastable state when F^ <T 
is a more delicate issue and requires the comparison of 
the spatial/temporal scale, F^/^, for critical thermal fluc- 
tuations to the typical expansion rate, H , expected in 



high-energy collisions. The small values of F^/*/i/ found 
in our calculations suggest that the phase transition is 
likely to proceed through spinodal decomposition rather 
than bubble nucleation for some fraction of all events (or 
some fraction of all rapidity bins). Our results certainly 
indicate that F^/"* is not materially larger than H. This 
suggests the familiar idealized Maxwell construction of 
equilibrium thermodynamics is not appropriate for the 
description of phase transitions in high-energy heavy ion 
collisions. In fact, the relevance of nucleation theory for a 
possibly first-order chiral phase transition in high-energy 
collisions appears questionable. The conditions F^ < T 
and If < Re leave at best only a small window of tem- 
peratures where nucleation might occur. 

Although current knowledge is insufficient to identify 
the order of the QCD phase transition in hadronic col- 
lisions with certainty, it is still important to search for 
experimental indications of supercooling. In addition to 
more familiar indicators, we have shown that the corre- 
lation length in the supercooled metastable state can be 
substantially larger than that expected in an equilibrium 
transition if the degree of supercooling is large. Whether 
this can lead to observable consequences depends sensi- 
tively on the expansion rate of the system, H. Small 
H will not support deep supercooling, and the relevant 
correlation length will remain close to h, which is small. 
On the other hand, large H will set causality limits on 
the size of the correlation length at the false minimum, 
which cannot exceed the horizon 1/ H. It remains to be 
seen if there is a "window" of energy or projectile/target 
combination for which the predicted increase in the corre- 
lation length can be observed. This speculation requires 
more detailed numerical studies of hydrodynamic expan- 
sion coupled to dynamical bubble nucleation (with real- 
istic initial conditions) before definite predictions can be 
made. 

The experimental observation of supercooling effects 
and spinodal decomposition would also be important as 
a matter of principle. Ideally, signatures of phase transi- 
tions should be order parameter related and should reveal 
properties of the equilibrium phase diagram. This is not 
the case for many of the signatures proposed for heavy 
ion collisions such as J/\l/ suppression and strangeness 
production. The interpretation of these signatures nec- 
essarily relies on theoretical estimates of reaction rates, 
and their connection to the thermodynamics of QCD is 
often far from clear. While the spinodal instability is not 
part of the equilibrium phase diagram, it is very close. 
Familiar experiments in condensed matter physics on a 
wealth of hysteresis phenomena (i.e. the analogue of su- 
perheating and supercooling) make it clear that spinodal 
instabilities can be studied on "macroscopic" time scales. 
The spinodal instability, if found, would probably repre- 
sent the most direct information which relativistic heavy 
ion collisions can provide regarding the QCD phase tran- 
sition. 

Not surprisingly, our results indicate that the time- 
scales for nucleation and expansion are not as well sep- 



arated in high-energy hadronic/heavy-ion colhsions as 
they are, for example, in the cosmological QCD phase 
transition. The assumption of a "well-mixed" phase, 
commonly employed in hydrodynamic calculations |2q ] 
and the opposite scenario of perfect capture of the order 
parameter in the metastable minimum [^ with subse- 
quent spinodal decomposition, must be regarded as rough 
qualitative pictures. Furthermore, the possibility of sig- 
nificant supercooling indicated by our results suggests the 
necessity of including the effects of non-spherical bubble 
configurations as well. Bubble growth and bubble perco- 
lation in particular are likely to play an important role 
in completing the transition. If the system evolves by 
strong supercooling followed by a rapid phase transition 
(due either to spinodal decomposition or nucleation of 
many small nonspherical bubbles which percolate), the 
constituent quarks will experience strong reheating due 
to the fact that they become essentially nonrelativistic. 
It would be interesting to seek experimental signatures 
for such "reheating" , which must be present if the quarks 
acquire mass instantaneously. 

In our view, a realistic and quantitative discussion of 
those points requires numerical studies of the dynami- 
cal evolution of quarks (which constitute the heat bath) 
coupled to the chiral field. While initial investigations 
along these lines have been performed, e.g. Refs. P,po|, 
the effects of dynamical bubble nucleation have not yet 
been considered. Moreover, the use of oversimplified ge- 
ometries did not allow for realistic velocity and temper- 
ature/density gradients, non-spherical bubble configura- 
tions, or bubble percolation. Work on these questions is 
in progress and will be reported elsewhere. 
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